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The Square Kilometre Array (SKA) will be operating at the time when several new large optical, 
X-ray and Gamma-ray facilities are expected to be working. To make SKA both competitive and 
complementary to these large facilities, thorough design studies are needed, focused in particular 
on imaging performance of the array. One of the crucial aspects of such studies is the choice of 
the array configuration, which affects substantially the resolution, rms noise, sidelobe level and 
dynamic range achievable with the SKA. We present here a quantitative assessment of the effect 
of the array configuration on imaging performance of the SKA, introducing the spatial dynamic 
range (SDR) and a measure of incompleteness of the Fourier domain coverage (Au/u) as prime 
figures of merit. 
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1. Introduction 

The demanding breadth of the science case and technical specifications of the Square Kilome- 
tre Array (SKA) makes the design of the array a complex, multi-dimensional undertaking (Jones 
2003). Although a number of factors may limit the performance of the SKA, including radio fre- 
quency interference, atmospheric and ionospheric effects (Perley 1999, Thompson et al. 1986), the 
array configuration is one of the most critical aspects of the instrument design (ISPO 2006), as it 
will be very difficult to modify or modernise the station locations after the construction phase has 
finished. 

Several approaches have been exercised in order to obtain an optimised array configuration 
(Conway 1998, Bregman 2000, Conway 2000a, 2000b, Kogan 2000a, 2000b, Noordam 2001, 
Cohanim et al. 2004, Bregman 2005, Kogan & Cohen 2005, Lonsdale 2005, Morita & Holdaway 
2005). Judgement of quality and fidelity of synthesised images is typically made by estimating 
the "dynamic range", the ratio between the peak brightness on the image and the r.m.s. noise in a 
region free of emission of the image (Perley 1999). High dynamic range is necessary for imaging a 
high-contrast feature, which is a key requirement for the SKA. SKA design goal for the synthesised 
image dynamic range is 10 6 at 1.4 GHz (Wright 2002, 2004). 

Alongside the dynamic range, one of the basic figures of merit (FoM) characterising imaging 
performance of an interferometer is the spatial dynamic range (SDR) quantifying the range of 
spatial scales that can be reconstructed from interferometer data (Lobanov 2003). The SDR of an 
interferometer depends on a number of instrument parameters, including the integration time of the 
correlator, channel bandwidth, and the coverage of the Fourier domain (wv-coverage). The latter 
factor poses most stringent constraints on the design, particularly for arrays with a relatively small 
number of elements. The quality of the wv-coverage can be expressed by the wv-gap parameter or 
Au/u, characterising the relative size of "holes" in the Fourier plane. Basic analytical estimates 
indicate that the SKA should have Au/u < 0.2 (Lobanov 2003) over the entire range of baselines 
to provide sufficient imaging capabilities and warrant that the SDR of the SKA would not be uv- 
coverage limited. 

The SKA configuration must provide nv-coverages that satisfy several key requirements de- 
rived from the prime science goals of the instrument: (i) good snapshot and deep imaging over 1 
degree field of view, (ii) high brightness temperature (Tb) for extended objects, (iii) dense core for 
transients/pulsar/SETI, and (iv) long baselines for milli-arcsecond imaging. The combination of 
these requirements with the benchmark figures for the dynamic range of continuum (10 6 ) and spec- 
tral line (10 5 ) observations poses a substantial challenge for the array design and for the antenna 
distribution in particular. 

A commonly used approach to designing the antenna configuration for an interferometric array 
relies on optimising the wv-coverage by minimising sidelobes or providing a desired beam shape 
(cf., Cornwell 1986, Kogan 2000a, 2000b). This approach assumes implicitly that the field of view 
is not crowded and the target objects are marginally resolved (so that the structural information 
can be recovered efficiently even if a substantial fraction of spatial frequencies is undersampled). 
Neither of these two assumptions will be correct for the SKA operating in the jiiJy regime. This 
complication requires additional constraints and considerations to be employed in order to warrant 
successful imaging of all spatial scales sampled by the interferometer. 
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This issue can be addressed effectively by ensuring a constant Au/u over the entire range of 
baselines, which will provide equal sensitivity to all spatial scales sampled by an interferometer 
and realize the full spatial dynamic range of the instrument. 

2. Spatial dynamic range of an interferometer 

Spatial dynamic range of an interferometer is determined by several major factors related to 
details of signal processing, design of the primary receiving element, and distribution of the collect- 
ing area of an array (Lobanov 2003). The maximum achievable spatial dynamic range (SDRp vX 
is given by the ratio of the field of view and the synthesised beam (HPBW). For an array composed 
of parabolic antennas 

SDR FoV «0.80|^, 

where, B max is the longest baseline in the array, T] a is the aperture efficiency, and D is the diameter 
of the antenna. However, SDRp v can be typically achieved only at the shortest baselines. The 
SDR is significantly reduced at the highest instrumental resolution, due to finite bandwidth and 
integration time, and incomplete sampling of the Fourier domain. Applying an averaging time of 
T a to interferometric data limits the maximum size a , of structure detected at full sensitivity to 
$a = c (v b s G)eT a B max ) _1 (Bridle & Schwab 1999), where v Q b s is the observing frequency and co e is 
the angular rotation speed of the Earth. Relating this size to the HPBW of the interferometer yields 
the maximum SDR that can be achieved at a given integration time: 

SDR T « (WeTa)" 1 « 1.13 X lO 4 ?" 1 . 

Bandwidth smearing due to a finite fractional bandwidth, Av, leads to a reduction in the peak 
response R& v = (1 + Av v /^hpbw) -1 ^ 2 (assuming a Gaussian bandpass and circular Gaussian 
tapering; see Bridle & Schwab 1999). The resulting limit on spatial dynamic range, SDRa v , can 
be approximated by the ratio v /#hpbw> which gives 

SDRav^Av-^av- 1 ) 172 ' 

with Rav < 0.75 typically assumed. Finally, incomplete sampling of the Fourier space yields 
SDRak < SDRpov- The magnitude of the SDR reduction can be expressed in terms of the "uv- 
gap" parameter, Au/u, that can be defined as follows: Au/u = (u2 — u\)/u\, where u\, U2 («2 > u\) 
are the wv-radii of two adjacent baselines (sampling respective structural scales 9\ j = l/wl,2, with 
6\ > 62)- Applying the assumptions used for deriving SDR Av gives a synthesised beam that can 
be well-approximated by a two-dimensional Gaussian. Then, for structures partially resolved at 
U2, the smallest resolvable size (or variation of the size) can be estimated by requiring that a dif- 
ference in visibility amplitudes V(u\) and V(u2) can be detected at a given SNR. This approach 
is similar to the one applied to determining resolution limits of an interferometer (c.f., Lobanov 
et al. 2001, Lobanov 2005), and for SNR > 1 it yields 2 M = (w/4)[ln2 ln(SNR)]- 1 / 2 . The 
ratio 62 /0\ can be represented by the term I + Au/u. With this term, the expression SNRa„ = 
exp[7T 2 (l+A M / M ) 2 /(161n2)] gives the relation between Au/u and the SNR required for detecting 
emission on spatial scales corresponding to the (u\ , u 2 ) range. In case of a filled aperture, for which 
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Ak/m -» 0, the corresponding SNR Au = SNR Ah / h=0 = exp[7T 2 /(16 ln2)] and SDR Al , = SDR FO v- For 
partially rilled apertures, the ratio SDRa m /SDRfov can be estimated from the ratio SNRo/SNRa„, 
which then gives 



Strictly speaking, the wv-gap parameter is a function of the location (u, 6) in the wv-plane (with 
6 describing the position angle), and it should be represented by a density field in the wv-plane. 
We will apply this description for making assessments of realistic uv-coverages obtained from our 
simulations. 

The 1/exp factor in the expression for SDRuv, calculated for two wv-points u\ and 112, es- 
sentially provides an estimate of a fraction of power that can be recovered between the respective 
angular scales (62, B\) from the sky brightness distribution. When an average value of Au/u over 
the entire wv-coverage is determined, the 1 / exp factor can be taken as a measure of ratio between 
the largest detectable structure and the primary beam (FOV) of individual array elements (under 
condition that the largest detectable size obtained from Au/u is smaller that the largest angular 
scale given by 1/ u m i n ). 

For an idealised, circular wv-coverage obtained with a regular array (i.e. logarithmic-spiral) 
with N stations organised in M arms extending over a range of baselines (B m i n ,B max ), the wv-gap 
can be approximated by Au/u rs (fi m ax/#min)^ — 1> with £, = M/N for baselines between antennas 
on a single arm, and ^ = l/N for all baselines. For instruments with multi-frequency synthesis 
(MFS) implemented, Au/u should be substituted by Au/u — AVmf s , where, Av m f s is the fractional 
bandwidth over which the multi-frequency synthesis is being performed. It should be noted that 
MFS will be not as effective improving Au/u for sources near the equator. 

In real observations, the actual SDR is determined by the most conservative of the estimates 
provided above. The simulations presented below will focus on SDRa„, assuming implicitly that 
the instrument is designed so that SDRa h > SDRa v and SDRa m > SDR T . Under this assumption, we 
investigate the relation between different array configurations as described by Au /u and resulting 
image properties described by the rms noise, dynamic range, and structural sensitivity of the array. 

3. Relation between Hv-coverage and imaging capabilities of an interferometer 

In order to provide a quantitative measure of the effect of array configuration on imaging perfor- 
mance of an interferometer, we simulate a set of an idealised array configurations, each providing 



over the entire range of u and 6. It should be noted that these configurations are introduced solely 
for the purpose of analysing the dependence of SDR on the wv-coverage, and they are not intended 
to serve as prototype configurations for the SKA 

We generate the test array configurations by considering an equiangular, planar logarithmic 
spiral and projecting this spiral on Earth's surface and determine the locations of individual sta- 
tions (i.e., latitudes and longitudes) using World GEOD system 1984 (Heiskanen & Moritz 1967). 
We apply Halley's third-order formula (Fukushima 2006), a modification of Borkowski's method 




Au 

— = const 
u 
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Table 1: Telescope settings for the generation of a range of visibility datasets. 



Parameter 


Value 


Frequency 


L band (1.4 GHz) 


Antenna 


SEFD= 335 K 1 " 


Bandwidth 


3.2 MHz 


No. of channels 


1 


Direction (J2000) 


00:00:00 +90.00.00 


Elevation limit 


12deg 


Shadow limit 


0.001* 


Start_time (IAT) 


22/08/2007 / 06:00 


Stop_time (IAT) 


22/08/2007 / 18:00 



Note: f - similar to the system independent flux density (SEFD) of a VLA antenna; $ - shadow 
limit is set such as not to constrain the wv-coverages obtained. 

Table 2: Input group of source components used for the generation of simulated visibility datasets. 



Source size 


RA 


Dec 


Flux density 




(J2000) 


(Jy) 


0.1' x 0.1' 


03:00:00 


88.00.00 


8.0 


1.0' x 1.0' 


08:00:00 


88.30.00 


4.0 


1.2' x 0.4' 


16:00:00 


89.00.00 


1.2 


3.0' x 1.0' 


21:00:00 


88.00.00 


3.0 


12.0' x 4.0' 


06:00:00 


88.00.00 


12.0 


42.0' x 14.0' 


18:00:00 


88.00.00 


42.0 


120.0' x 40.0' 


00:00:00 


90.00.00 


120.0 



(Borkowski 1989), to transform Cartesian to geodetic coordinates. This method is sufficiently pre- 
cise; the maximum error of the latitude and the relative height is less than 6 /xarcsec for the range 
of height, —10 km < height < 30,000 km, and is stable in the sense that it converges for all co- 
ordinates including the near-geocentre region and near-polar axis region. We then transform these 
coordinates to obtain ECEF (Earth Centered Earth Fixed) coordinates using a glish script, and pro- 
duce an input array configuration file. The actual dimensions of the simulated array are chosen 
such that baseline non-coplanarity is negligible and no w-term correction is required. This does not 
affect conclusions obtained from the simulations, which are generic and can be applied to evaluate 
imaging performance of any interferometric array of arbitrary configuration and extent. 
In order to generate wv-coverages satisfying the condition 

Au . 

— (u, 0) = const, 

we place the array center at the North Pole, and consider a fiducial field centered at a 90 deg declina- 
tion (which yields circular wv-coverages that can be described by a single value of Au/u). The corre- 
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sponding visibility datasets have been generated in aips++ (ver 1.9, build 1556; casa . nrao . edu) 
Imaging has been done with MPS package (www . aips . nrao . edu). Several Fortran programs 
and glish scripts have been developed and used to automate the process and provide a structured, 
uniform, repeatable and robust processing. 

We generate datasets for full track observations (IAT 6:00-18:00 hrs), using a range of inte- 
gration times at a frequency of 1.4 GHz and with a 3.2 MHz bandwidth (see Table [[]). Antenna 
sensitivity (listed in Table [I]) are assumed to be similar to a dish of Very Large Array (VLA). The 
simulator does not take into account the effect of the primary beam on imaging. An elevation limit 
(horizon mask) of 12 deg and a shadow limit (maximum fraction of geometrically shadowed area 
before flagging occurs) of 0. 1 per cent was introduced. 

We produce a number of array configurations corresponding to different values of Au / u and 
use them to obtain visibility datasets for the chosen model brightness distribution (see Table ||). We 
repeat this exercise for all generated array configurations and perform identical pipeline analysis, 
to ensure a self-consistent comparison of basic characteristics of the resulting dirty and CLEAN-ed 
images. 

3.1 Simulations 

This section describes simulations and analysis of array configurations providing equal wv-gap 
at all baseline lengths. 

We have tried two different algorithms for generating array configurations with a given Au/u = 
const for all baselines. Both algorithms have employed a logarithmic spiral geometry, but differed 
in the realisation of changing the characteristic value of Au/u from one configuration to another. 
In the first alogrithm, the total number of stations was kept constant, while the baseline spread was 
gradually increased. In the second approach, the baseline spread was kept constant, and changing 
Au I u was achieved by changing the total number of antennas in the array. We have found that the 
second approach is superior for maintaining a constant noise level for different array configurations, 
and we adopted it as the basis for our simulations. In this set of simulations, the largest baseline 
length, fi max = 5 km is kept constant and different values of Au/u are realized by varying the total 
number of antennas N located on a single arm of an equiangular, logarithmic spiral. Fig. [l] gives an 
example of an simulated snapshot wv-coverage generated using this approach. 

We consider only baselines to the station at the origin of the array and generate a range of 
visibility datasets to probe the wv-gap parameter from 0.19 (N = 50) to 0.01 (N = 640). The simu- 
lated intergation time is set to 1 second, and the resulting sampling times range from 12.5 seconds 
(one wv-point every 12.5 seconds) for N = 50 to 160 seconds for N = 640. With these settings, we 
produce visibility datasets for full track wv-coverages for each of the array configuration. 

>From the simulated data, we produce both dirty and CLEAN images, of 8192 pixels x 8192 
pixels in size, with a pixel size of 3 arcsec. The dirty images obtained for several different simulated 
configurations are shown in Fig. [2[ The respective interferometer beam cuts are shown in Fig. ||. 
The CLEAN images were obtained using AIPS task IMAGR applied to 499 facets spread across 
a ~ 50 square degree field. It should be noted that the CLEAN algorithm has been applied non- 
interactively (i.e., with the same set of CLEAN parameters for all datasets and without varying 
them during the CLEAN-ing), and only a fraction of the total flux density has been recovered 
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Figure 1: An example of a simulated 2-minute long snapshot Mv-coverage outlining relative locations of 
antenna stations. The largest baseline length, B max = 5 km, and the number of antennas is N = 50; In the 
simulations, only baselines to the central station are considered, resulting in Au/u = 0.19 for a full track 
observation. 



after 200,000 iterations, which points out to likely limitations due to application of gridding and 
deconvolution. 

The resulting facet images were stitched together using AIPS task FLATN to create a single 
final image. We used uniform weighting and the 3D option for the w-term correction throughout 
our analysis (Cornwell & Perley 1992). In order to assess the effect of deconvolution and gridding 
of wv-data on the results of the simulations, both the CLEAN and the dirty images were used for 
estimating several basic FoMs. We use AIPS task IMSTAT to determine the r.m.s. noise levels in 
each case. The results of these estimates are compared in Fig. ||. 

One can see that the dependence of the rms noise (and reciprocally, the dynamic range) on the 
wv-coverage changes at Au/u 0.03. At smaller values of Au/u, the wv-coverage does not have a 
strong effect on the flux density recovered from the visibility data (this holds true for both dirty and 
CLEAN images). The requirement of 

Au/u < 0.03 

for the entire range of baselines can therefore be used as a benchmark requirement for designing 
the SKA configurations that would minimise the effect of Mv-coverage on reconstructing the sky 
brightness distribution. It should be noted that this conclusion provides a strong benchmark, im- 
plying implementing wv-coverages with wv-gaps larger than 0.03 will limit the imaging capabilities 
of an array even with although the present day convolution and imaging algorithms. 

Further, more detailed investigations may be required in order to refine this conclusion and 
assess the full range of effects that may potentially affect the dynamic range and structural sen- 
sitivity derived from images obtained with different wv-coverages and different values of Au/u. 
One possibility would be to use more elaborate and more realistis sky models coming from source 
simulations (e.g., O'Sullivan et al. 2009). 
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Au'U- 0.M5 



Au u- 0X112 




r 



AU u 1.1 Of* 



Figure 2: Dirty images obtained from the simulated data. Each panel is marked with the value of Au ju for 
the respective simulated dataset. 



4. Practical issues 

We illustrate now a practical way to use the Au/u parameter for evaluating an arbitrary array 
configuration. The evaluation is based on analysis of a Mv-coverage arising from a specific obser- 
vation or a set of observations (for instance, a number of snapshots on targets covering a specified 
range of declinations and/or hour angles). 

Since almost all real wv-coverages deviate from circular symmetry, the Au/u FoM should be 
represented either by a two-dimensional distribution or by an average of that. The averaging can 
be made azimuthally (providing a profile of Au/u as a function of wv-distance) or both azimuthally 
and radially (giving a single value description of a given wv-coverage). For the averaged quantities, 
the respective dispersions can be used to quantify inhomogeneities in the i<v-coverages. 

For the purpose of a better graphical representation of the wv-gap distribution, a definition 



8 



SKADS: Array Configuration Studies 



,au/u = 0.190 




2000 4000 6000 





Au/u = 0.045 
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2000 4000 6000 



Au/u = 0.022 




A u/u = 0.092 




Au/u = 0.060 



6000 8000 



2000 40 
pixel n 


600 D BO 


DO 


2000 40 
pixel n 


6000 80 




Au/u = 0.036 


0.06 
0.04 




Au/u = 0.030 






E 

$ 0.02 
0.00 
-0.02 
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iiiimiiiiMi'ilkliiii'iiiiill 


2000 40 
pixel n 


6000 80 


00 


2000 40 
pixel n 


6000 80 




Au/u = 0.019 


o.oe 

0.04 




Au/u = 0.014 


lillllllilill | iBfcn 

2000 40 


ihMi^lllLiliiiilillii'i) 
lipipiipiflllipiip 

6000 80 


E 

% 0.02 
0.00 
-0.02 

DO 


2000 40 





Figure 3: One-dimensional profiles of restoring beams obtained from selected datasets. Each panel shows 
the value of Au/u corresponding to the respective simulated dataset. The profiles are truncated at 0.06 of the 
beam peak in order to emphasise the difference in the sidelobes resulting from different wv-coverages. In 
each panel, the red curve represents the 30-pixel boxcar smoothing of the respective beam profile. 



Au/u = (u2 — u\)/u2 («2 > u i) can be applied. This definition changes the range of Au/u from 
[0,°°] to [0,1], with corresponding to Au/u from two identical baselines. The maximum value 
Au/u = 1 is realized everywhere outside the area covered by the observation (for which U2 = °° can 
be assumed). 

A feasible approach to determine the figures of merit based on Au/u is as follows: 

1. Grid the wv-data into N sectors, where the width of each sector is determined by the observing 
scan-length (the duration of a snapshot observation, or a typical length of a single scan in a 
synthesis observation). 

2. For each individual sector (described by its central position angle, 0, and width 2n/N and 
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Figure 4: Relative changes of different FoMs obtained from the dirty (left anel) and CLEAN (right anel) 
maps. The normalization is done with respect to the FoM values determined from the dirty image with 
Au/u — 0.01, by normalizing its beam area to unity, rms noise to the value of 2, peak surface brightness to 
the value of 3, and dynamic range to the value of 6. The peak surface brightness recovered from different 
simulation runs remains nearly constant. The rms noise first decreases rapidly for smaller values of Au / u, but 
the rate of this decrease is then remarkably reduced, implying that the wv-coverage does not affect strongly 
the image noise and dynamic range for array configurations with Au/u < 0.03. Note that the dynamic 
range is determined for a feature with the size substantially smaller than the largest angular scale for the 
wv-coverage with Au/u = 0.19, and thus the calculations are not affected by variable minimum Mv-spacing 
in different simulation runs. 

containing M wv-points), determine Au/u(u, </>;) for all M wv-points. As was pointed above, 
Au/u(u, fa) < 1 within the range u\ <u <um and it is equal to unity elsewhere. 

3. Plot the combined results for all sectors in polar or rectangular coordinates or perform aver- 
aging in azimuth (^f (m,^))^ and/or azimuth and radius (^f («, <P))u,<j>- 

A detailed description of the calculaton algorithm is given in Appendix 1. Appendix 2 sum- 
marizes a proposed set of FoM that includes the wv-gap parameter and can be used for evaluation 
and optimization of the SKA configuration. 

4.1 Examples 

In order to represent a two-dimensional distribution Au / u, the Voronoi tessellation can be ap- 
plied to combine the Au/u calculated for the individual sectors. The Voronoi tessellation ("Voronoi 
diagram") uses partitioning of a plane with n points into convex polygons such that each polygon 
contains exactly one generating point and every point in a given polygon is closer to its generating 
point than to any other (Okabe et al. 2000). The Voronoi diagram is sometimes also known as a 
Dirichlet tessellation. The cells are called Dirichlet regions, Thiessen polytopes, or Voronoi poly- 
gons. We use the algorithm by Okabe et al. (2000) built in MATLAB and GNU Octave (version 
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Table 3: Table showing values of Au/u from azimuthal-radial averaging and the corresponding dispersions 
for all panels with numbers 1 through 6 in Figs. 6 and 7. 



Example wv-coverage 


(Au/u) 




(1) 


Simulated log-spiral: core baselines 


0.824 


0.071 


(2) 


Simulated log-spiral: all baselines 


0.219 


0.244 


(3) 


Simulated skipped spiral^: all baselines 


0.203 


0.246 


(4) 


VLA snapshot 


0.072 


0.084 


(5) 


GMRT short observation 


0.016 


0.043 


(6) 


GMRT nearly full 12 hr synthesis 


0.012 


0.035 



Note: t - a multi-arm, logartihmic spiral configuration, with concurrent antennas removed from all but one 
of the arms of the spiral, in order to reduce the number of redundant baselines as much as possible. 

3.2) to decompose the discrete set of — (uv, 0; (i = l,N) values into density plots that can be used 
as a graphical representation of the wv-gap distribution. 

We exemplify this approach by considering six different wv-coverages from simulated and real 
observations (Fig.|5]). The resulting two-dimensional distributions of Au/u are shown in Fig.|6[ One 
can either analyze these distributions directly, or produce azimuthally-averaged radial profiles, or 
describe these Mv-coverages by a single FOM obtained from double, azimthal-radial averaging. 
Table |3| shows values of Au /u from azimuthal-radial averaging and the corresponding dispersions 
for all panels with numbers 1 through 6 shown in Fig. 0. 

The radial profiles obtained by azimuthal averaging of the density plots are shown in Fig. 7. 
One can see immediately, from the profiles in Fig. 7 and the dispersions of the mean in Table ||, 
that Au/u varies strongly with the baseline length. This dependence should be minimised during 
the design of the SKA configuration, taking into account the geographical location of the array, the 
different types of observations (e.g., snapshots, full track synthesis and so on), and the range of 
declinations for which the imaging performance of the SKA should be optimised. 
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Figure 5: Six different Mv-coverages used to exemplify determination of the Au/u figures of merit: (1) a 
wv-coverage from a simulated log-spiral, baselines to the core station included; (2) the same simulation, 
all baselines included; (3) a w-coverage from a simulated skipped-spiral configuration (aimed to provide 
a nearly constant Au/u and to satisfy the collecting area distribution required for the SKA); (4) a VLA 
snapshot Mv-coverage; (5) a GMRT Mv-coverage from a short observation; (6) a GMRT Mv-coverage from a 
nearly complete, 12-hour synthesis. 
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5 10 15 -15 -10 -5 

u (kX) u (Kk) 



Figure 6: Density plots of -j-(u,^>) distribution calculated from the Mv-coverages shown in Fig 5. Note that 
for very sparce Mv-coverages, the tesselation algorithm may introduce artefacts. This can be illustrated by the 
apparent asymmetry of Au ju distribution seen in panel 1 produced for a symmetric Mv-coverage with very 
few Mv-points at large t/v-radii). More refined approaches to calculating the density fields from Mv-coverages 
may be needed for such cases. 
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(1) ' Simulated log-spiral: core baselines' 



10 15 20 25 

radial-distance (kilo-lambda) 



{2} 'Simulated log-spiral: all baselines' 
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(3) 'Simulated skipped log-spiral: all baselines' 



10 15 20 25 

radial-distance (kilo-lambda) 



(4) 'VLA snapshot' 



50 100 150 200 250 300 
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(5) 'GMRT short observation' 
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(6) 'GMRT (nearly) full 12 hr synthesis' 



2 4 6 8 10 12 14 16 18 20 
radial-distance (kilo-lambda) 



Figure 7: Azimuthally averaged profiles of Au/u from the density plots of ^(w, 0) shown in Fig. 7 for all 



six cases. 
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Appendix 1: An algorithm for calculating Au/u from Hv-coverages 
Initial settings 

Let us suppose we have a file containing visibilities or a wv-coverage corresponding to an obser- 
vation with a target array configuration. The wv-coverage is given in polar wv-coordinates (u, 6) 
and spans a range (« m in,"max) of wv-distances. Note that the wv-coverage can be completely arbi- 
trary: i.e., it can have asymmetries and gaps along certain position angles — this does not affect 
the calculation. 

The calculation must be done for N individual intervals in position angle. The number N is 
determined by the dynamic range specification for the array so that N = SAT? 1 / 2 . This corresponds 
to a width of Ad = 180/SNR 1 / 2 degrees for an individual interval and, respectively for a time 
integration of At = 720/SNR 1 / 2 minutes. For a target SNR ->■ 10 6 specified for the SKA, one can 
assume (for convenience of calculation) 

N = 720, At = lmin, Ad =0.25deg. 

The wv-tracks on individual baselines should either be generated with an integration time At or 
averaged to this time (see Step 2 of the calculation). 

For convenience of notation, let us denote 

Su = Au/u 

in the following description of the calculation algorithm. 
Calculation of Au/u 

Step 1. Break the Fourier plane (u, 6) into N intervals in B (0-intervals) and N intervals in u (u- 
intervals), thus forming mNxN grid in (u, 6) with individual grid cells (w-, with i = \,N 
and j = l,N. Note that 6 € (0°, 180°) and u G (0,w max ). 

Step 2. If needed, average the wv-data to At minutes to ensure having one wv-point per baseline per 
(9-interval. 

Step 3. In a given ^-interval, 6j (j = l,N), add the zero-spacing point (w = 0) to the wv-points 
falling into this interval and calculate (ungridded) values of 8u from 

8uk = (uk-Uk-i)/uk, 

where k = l,M and M is the total number of wv-points falling within this d -interval. 

Step 4. For the given ^-interval, dj, map the values 8uk to individual grid cells (w-, d'j) as follows: 

VmJ 6 (uk-\,Uk) : 8ui = 8uk, k=l,M,i=l,N. 

• For each grid cell, check the following two conditions: 

— > If an j-th grid cell contains one of the ungridded values , 8uk, then set its wv-gap value to 
8uk — that is, no interpolation between the 8ut and 8uk+\ is needed for this grid cell. 



17 



SKADS: Array Configuration Studies 



— > If an j-th grid cell contains several ungridded values of Su, set its Su value by averaging 
these 8u values. 

• After mapping the last (M-th) ungridded value of Su, check the following condition: 

— >■ If um < "max, assign 8u = 1 to the grid cells falling within the interval (uM,u max ), starting 
from the next grid cell after the cell containing the wv-point um- 

As a result of this Step, the given ^-interval is divided into N grid cells in u each assigned a 
value of 8u. 

Step 5. Repeat Steps 3 and 4 for all ^-intervals and populate the entire (V, 6') grid with wv-gap 
values Suij, where i = l,N is the index in w-intervals and j=l,N is the index in 0-intervals. 

Step 6. Plot the gridded distribution 5w,-j for visual inspection, in the rectangular coordinate system 
(u, G). No smoothing in u or 6 is required. 



Step 7. Calculate N azimuthal averages of 8u 



1 N 

(Hi = 77 L Sltij 
ly 7=1 



and their variances 

N 



of = \y{8u i , j -{8u) i ) 2 . 



Step 8. Plot the resulting one-dimensional distribution (u' h (Su)^, with the respective dispersions, <J; 
as errorbars. This gives a radial profile of the wv-gap. 

Step 9. Calculate the integral FoM value of 8u and its variance from 

1 N 

(Au/u) u , e = -Y(Su)i, 
1 N 

2 _ V 2 
G Au/u — 77 2- °i ■ 
n (=1 

Various refinements to the procedure described in above can be made. One possibility, for instance, 
is to calculate Au/u within an angular sector with the width of Ad, then shift this sector by an angle 
corresponding to the observation intergration time and repeat this procedure until the entire range of 
the position angles is covered. The resulting distribution of Au/u can be then averaged or smoothed 
over Ad. 
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Appendix 2: A set of figures of merit for SKA configuration evaluation 

This appendix describes a proposal for an initial set of figures of merit (FoM) for evaluating imag- 
ing performance of different array configurations for the SKA. The set of FoM proposed includes 
general performance metrics and does not include specific requirements on array configuration 
coming from individual SKA Key Science Projects. The critical issues for evaluating imaging 
performance provided by a given array configuration are: 

1. Providing optimal shape of the point-spread-function (PSF). 

2. Minimizing sidelobes of the PSF. 

3. Evaluating sensitivity to all spatial scales sampled by the SKA. 

Table Al sumarizes a basic set of seven FoM suitable for evaluating the three conditions listed 
above. The FoM listed in Table 1 are decribed in polar coordinates in the wv-plane (u,6) and image 
plane (r,d). 



Table Al: A set of figures of merit for SKA configuration evaluation 



FoM 


Description 


Name 


Goal 


Point Spread 


Major Axis 


BMA, 






Function (PSF) 


Minor Axis 


BMI, 




^min/^maj ^ 1 




PA. of Major Axis 


PAN, 


^maj 






PSF Shape 


PSF, 


B(r,9) 


B(r, 6) — > Gaussian 


PSF Sidelobes 


Maximum Positive Sidelobe 


MPS, 


o+ 


a+ ^0 




Maximum Negative Sidelobe 


MNS, 


o 


a ^0 




Sidelobe RMS 


RMS, 


Orms 


Orms -> 


UV Gap 


Integrated Value 


UVG, 


Au/u 


(Au/u) Uj e ->0 




Dispersion 


UVD, 


Au/u(u, 6) 
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